Testing dark matter warmness and quantity via the reduced relativistic gas model 
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We use the framework of a recently proposed model of reduced relativistic gas (RRG) to obtain 
the bounds for O's of Dark Matter and Dark Energy (in the present case, a cosmological constant), 
taking into consideration an arbitrary warmness of Dark Matter. An equivalent equation of state has 
been used by Sakharov to predict the oscillations in the matter power spectrum. Two kind of tests 
are accounted for in what follows, namely the ones coming from the dynamics of the conformal factor 
of the homogeneous and isotropic metric and also the ones based on linear cosmic perturbations. 
| The RRG model demonstrated its high effectiveness, permitting to explore a large volume in the 

space of mentioned parameters in a rather economic way. Taking together the results of such tests 
as Supernova type la (Union2 sample), H(z), CMB (R factor), BAO and LSS (2dfGRS data), we 
confirm that ACDM is the most favored model. At the same time, for the 2dfGRS data alone we 
found that an alternative model with a very small quantity of a Dark Matter is also viable. This 
output is potentially relevant in view of the fact that the LSS is the only test which can not be 
affected by the possible quantum contributions to the low-energy gravitational action. 

PACS numbers: 98.80.-k, 98.80.Cq, 98.80.Bp, 98.80.Es 
^ ^ ■ Keywords: Warm Dark Matter; Ideal Relativistic Gas; Cosmic Perturbations; Quantum effects of vacuum. 
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P 1 I. INTRODUCTION 

According to the standard estimates of the energy balance of the present-day Universe, the relative energy densities 
of the Dark Energy and Dark Matter (DM) are close to =0.7 and = 0.25, respectively, while the visible 

{SJ ' (more precisely, baryonic) matter is represented by a modest less than 5% of the total energy density [l], 0. Both 
^ . Dark Energy and DM represent somehow mysterious components, but for different reasons. The main candidate to 
be Dark Energy is a cosmological constant (CC), however the cosmological constant problem is considered as one of 
the most difficult problems of the modern fundamental physics. The main CC problem is due to the conflict between 
the overall value of CC which is likely causing an acceleration of the present-day universe, and several much greater 
contributions to it which can be evaluated at the particle physics scale (see Q for a review). At the moment, no 
solution of this problem is known. On the other hand, there is another difficulty, which is known as a problem of 
coincidence. The problem is to explain why the universe started to accelerate only recently. This problem can be 
solved or at least alleviated in the two distinct ways. First, one can replace the CC term by some artificially designed 
substance (e.g. quintessence) which can adjust its energy density according to the expansion of the universe. Indeed, 
this approach is essentially worsening the situation with the main CC problem, because the amount of the requested 
fine-tuning becomes even greater and goes beyond the original one, which is due to the hierarchy between the particle 
physics scale and cosmic scale Q. The second approach is much less explored, despite it looks much simpler. It 
C$ ■ assumes changing the energy balance of the universe such that in the new model becomes closer to unity. In this 
case the acceleration starts earlier and the coincidence problem becomes partially resolved. 

The unique way to increase Q° is by the expense of Q.% M , and from the first sight this idea does not look realistic. 
The well-known reason is that the DM is required by a well-established set of observational data (see, e.g., books 
and reviews At the same time it is interesting to verify the mentioned possibility without prejudice, especially 

taking the possible warmness of the DM into account [H, [lij]. The exploration of the wide set of observational 
data is a complicated task. In particular, one has to choose such a model which could take care about possible 
DM warmness and, at the same time, keep the amount of requested calculations under control. It is also desirable 
to be able to implement some new physical input into this model, for the case we will like to make it at the next 
stages of investigation. All in all, we need sufficiently simple, but yet physically reasonable model for a warm DM 
(WDM). Fortunately, a very useful simplified model of this sort, called reduced relativistic gas (RRG), has been 
suggested recently in [lj| [Tfl]. An important historical note has to be done. After the first version of this paper has 
been submitted, we have learned from [l8| that the equivalent simplified model, interpolating between radiation and 
matter epochs, has been used by A.D. Sakharov in the famous work (l7j | . where the oscillations in the matter power 
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spectrum were discovered for the first time. 

The idea of RRG is to treat WDM as an approximately Maxwell-distributed ideal gas of massive particles. As far 
as the large-scale description of the Universe concerns mainly the equation of state of the DM, it is not supposed to 
be sensitive to the non-collisional nature of the DM particles and to their anisotropic distribution at the astrophysical 
scale. Hence, there can not be much difference between the RRG description from one side and the one based on 
relativistic Fermi distribution from another one. The same is expected to be true, up to some extent, for the more 
complete approach based on the Boltzmann- Einstein system of equations @ . 

The expression "approximately Maxwell" concerning the RRG model means we take an approximation where all 
DM particles have equal kinetic energies. As it was shown in (l5j |. qualitatively and numerically the equation of state 
suffers only a very small modification from this approximation. At the same time, the model based on RRG admits 
simple analytical form for the equation of state and even enables one to analytically solve the Friedmann equation. 
The analysis of cosmic perturbations in the framework of RRG model and LSS data provided (l6j an upper bound 
on the DM warmness which is very close to the one obtained from much more complicated analysis based on the 
standard Boltzmann- Einstein system [L9U22I ] . Thus, the RRG model represents a really useful tool for exploring the 
new cosmological models with the WDM. At the same time, it looks interesting to check this model on such sets of 
observational data as CMB, BAO and others. 

There is also a special motivation for us to explore the possibility of an alternative values of Q-dm with the WDM, 
which is related to the potentially relevant quantum effects of vacuum. The backgrounds of the cosmological and 
astrophysical applications of the theory with renormalization-group based quantum corrections were established in 
[2~il . In the recent paper [27j it was shown, from a general considerations based on covariance and dimensional 
arguments, that the form of these quantum corrections can be defined up to a single free parameter v. The recent 
analysis of the galaxies rotation curves in [23j | in the theory without DM, but with the mentioned quantum corrections 
to the Newton law was surprisingly successful. It was shown that such quantum corrections can compete with CDM 
in explanation of the rotation curves. The small amount of sufficiently warm DM does not affect these results and, 
therefore, we can think about weakening the standard restrictions to the DM in the theory with quantum corrections. 
On the other side, the cosmological model based on the same renormalization-group quantum corrections has been 
considered in [24,25]. The analysis of cosmic perturbations in this model [25] (see also [26[) demonstrated that the 
quantum contributions almost do not affect the power spectrum, such that their impact on the LSS (2dfGRS) data is 
very weak. Therefore, if such quantum corrections do exist, the LSS data alone should manifest such possibility even 
in the zero-order approximation, that means without taking into account the quantum effects. At the same time, 
we should expect the quantum corrections to be relevant for other test, especially for lensing, Supernova, H{z) and 
CMB. In order to take this into account, one has to consider the 2dfGRS test separately. The positive result for the 
zero-order approximation is the case when the LSS part alone should be compatible with an alternative model, with 
smaller £Idm> while the full set of tests should uniformly indicate to the ACDM as unique possible option without 
taking quantum effects into account. This is the only one possible output which can actually leave the chance for 
relevant effects of quantum corrections at the astrophysical and cosmological scale. 

We will present in this paper the RRG model of a WDM, characterized by the dark matter fractional density 
measured today fidmo an d also by the parameter 6, which specifies the degree of warmness of the DM component. The 
model will be confronted against background observational data (Supernova type la, Baryonic Acoustic Oscillation, 
Cosmic Microwave background R parameter and the age of the universe as function of z, H(z)) and the matter power 
spectrum at the linear regime. We will see that the background tests, in general, favor a scenario of the ACDM model, 
with b rj 0, and Sldmo ~ 0.25, even if for a given specific tests some particularities may appear. However, the matter 
power spectrum alone favors a model with a very small value of fidmOi giving the room for some degree of warmness. 
In view of the discussion presented above this result can be seen as a sign for a possibility of a new concordance 
scenario if the quantum corrections are taken into account. 

The paper is organized as follows. In Sect. 2 we present a brief necessary review of the RRG model. The contents 
of this section is essentially the same as [16[ and we include it only for making all considerations in the paper more 
consistent. In Sect. 3, we discuss linear perturbations of matter including baryon component. The observational 
data which are used in our analysis are presented in Sect. 4. Finally, in Sect. 5 we draw our conclusions and discuss 
further perspectives for exploring an alternative concordance model with quantum corrections. 

II. BACKGROUND NOTIONS OF THE RRG MODEL FOR WDM 



In this section we present the necessary elements of the RRG model [l5j, ]1 
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A. Equation of state and zero-order cosmology 

Consider the cosmological model which is based on the relativistic gas of massive particles representing WDM, plus 
baryonic matter, radiation and cosmological constant. As a simplification, it is assumed that the relativistic gas is 
composed by particles with equal kinetic energy, or equal speed /3 = v/c. By elementary means one can obtain the 
following equation of state, which links the pressure P and energy density p of such gas, 
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In this formula e is the kinetic energy of the individual particle, e = mc 2 /\/l — /3 2 . Furthermore, pd = AqO- + z ) 3 
is the mass (static energy) density. One can use one or another form of the equation of state ([TJ , depending on the 
situation. For example, in order to explore the evolution of the conformal factor a = a(t) in this model, one can take 
the second expression in ([1]) and using the conservation law find that the relative energy density for the relativistic 
gas is given by the expression 



tW(a) = Jl + ^, (2) 

a 3 V 1 + b 2 V a 2 

where fidmo is the present-day dark matter density, and b is related to particles speed. Indeed, b rs means that the 
particles are nonrelativistic. Furthermore, for small velocities one can easily obtain the relation b « /3 (l6| . Finally, 
by using the equations presented above and the ones for other components, the Hubble parameter can be presented 
in the form 



H 2 ( a , n dm0 , n b0 , b, n r0 ) = Hi \ ndm ° J i + _1 + 0*° + + n Aa) ] , (3) 

LaVl + b 2 V a 2 a 4 J 

where we assume that 

fiAO = 1 — "dmO — Slfco — "rO 

and use the value Q r oh = 2.42 x 1CP 5 . In this work we consider only the space-time manifolds with the flat space 
section. 

If taken alone, the RRG model provides a natural interpolation between the radiation-dominated and matter- 
dominated universe regimes. Due to the adiabatic expansion, the relativistic gas is cooling down and the universe 
can look as filled by ultrarelativistic gas at the early stage of expansion and almost as a dust at the final stage of this 
expansion. One can find more details about the RRG model in [l5| and [l^ . 

In order to illustrate the effect of the warmness of the matter content of the universe, in Fig. [T] we displayed the 
deceleration parameter q(z), the equation of state w(z) and the age universe in Hq units. One can see that for b ~ 0.1 
the function plot goes away from the ACDM plot with (6 = 0). The state function is the function more sensitive to 
the value of parameter b. 



III. EQUATIONS FOR LINEAR PERTURBATIONS 



Following the approach developed in our previous paper [161 ] . the dynamics of density perturbations has been 
analyzed in the synchronous gauge. The calculation of these perturbations is rather a standard routine, the only 
exception is the variation of the equation of state ([TJ , where one has to use the first form of this equation, taking 
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The reason is that in the framework of RRG model one has to provide kinetic energies of all particles to be equal and, 
therefore, we have no right to change the ratio mc 2 /e. In other words, the variations of the energy density and of the 
rest energy density, 5pd and dp, should b e al ways proportional. By using this rule we arrive at the following system 
of equations for the linear perturbations lit 
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FIG. 1: The deceleration parameter (left), equation of state (center) and age universe (right) for different values of the parameter 
6 = 0, 0.01, 0.2 (from top to bottom of the first plot). Here w(z) is an accumulated parameter corresponding to the sum of 
the cosmological constant term and the matter content. In all cases here we used the values Qi,o — 0.04 and Q<im0 = 0.25. The 
ACDM model corresponds to the value b = 0. 



where v is the peculiar velocity, <5& is the baryonic matter density fluctuation and 8dm is the relative dark matter 
density fluctuation. Other relevant definitions are 



/ = /(«) = a "\j ^dm + -Jp + 1 — fibO — ^dmO , 

1 - r(ffl) l+b 2 

w = w(a) = and r{a) = - - , (6) 

where (Idm is given by the Eq. © . The system above differs from the CDM model, because in our model the peculiar 
velocity and the state function w(a) are non-zero. Additionally, the CDM model interacts only through gravity and 
can be considered as a pressureless perfect fluid, unlike our model has a component of pressure given by (jl]). The 
solution of the system (5) w as p erformed using a Mathematica-based code. To determine the power spectrum we use 
the BBKS transfer function [28| . It is important to note that we also use the transfer function described in references 
[2fjl I30I ] which are adapted for WDM particles, which depends on some specific inputs, obtaining the same results as 
when the BBKS transfer function is employed. 

For the scale invariant spectrum, favored by the primordial inflationary scenario, the BBKS transfer function is 
given by 

= ^ (1 + 2.34g) + + 16 + 5 + 6 7 4] -i/4 (7) 
2.34g 

where q(k) — — — - and V = fl m Q h cxp ( — flbo — 



hT Mpc 1 \ !] m o 

We have used the function T(k) to impose the initial conditions of the system (5). The power spectrum is defined 
as usual, at z — we have 

P(fc) =\ S (k)\^AT(k)\ 9 -^] 2 k, (8) 

where A is a normalization constant of the spectrum, which can be fixed from the spectrum of anisotropy of the 
cosmic microwave background radiation. The expression for g(f2) is given by 

= 2(Sl 4 /7 + i.oi(f + 1) -0.75) ^ 

The analysis of the linear density perturbations enables us to obtain the upper bound for the warmness of the DM 
in a very economic way [l6| . In Fig. 1 we compare the 2dFGRS data with the power spectrum of our model for the 
usual ACDM energy balance case. We used the best fit obtained using the \ 2 m the next section. One can see that 
if b « 10~ 3 , the predicted power spectrum is outside of the region of data shown in the plot. A detailed analysis of 
observational constraints is given in section 5. 
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FIG. 2: Power spectrum using n 60 = 0.044 and fidmo = 0.24, for three different values of b = 10~ 4 , 5.0 x 10 -3 , 10~ 2 (from 
bottom to top). It is easy to see that the last two values for b provide the linear power spectrum, at small scale, which is 
incompatible to the 2dFGRS data. 



IV. CONSTRAINTS FROM OBSERVATIONAL DATA 



Nowadays there is a large number of observational tests focused on different aspects of the cosmological models. 
It is expected that the crossing of these different observational tests may lead to an expressive constraint on the 
free parameters of a given model, in some cases confirming the viability of the proposed scenario or ruling it out. 
Here we consider essentially two kinds of tests: those focusing on the background functions, specially H{z) and also 
the perturbative tests, which involve more profoundly the description of the matter/energy content of the model. 
Concerning the tests related to the background behavior only, we will consider the supernova type la data (Union2 
sample), BAO, the position of the first acoustic peak in the CMB spectrum (the R parameter) and the age of the old 
objects. In what concerns the perturbative tests, we will restrict ourselves to the matter power spectrum associated 
to linear perturbations. 

The analysis of statistics begins with the x 2 functions constructed following the general expression 

N f n th n ob\2 

x 2 = J2 {pi ~J i] , (10) 

i—n 1 

where N is the total number of observational data, p* 71 are corresponding theoretical predictions and p° h represent 
the observational values, with an error bar given by a%. The probability distribution function is constructed from the 
X 2 function as 

P(x j ) = Ae^ 2 ^''/ 2 , (11) 



where x J are the set of free parameters of the model and A is a normalization factor. Within the RRG model, there 
are two free parameters, Od m o and b. In what follows we describe the results for each of these tests. 
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A. The CMB Shift Parameter 

We used the shift parameter, R, which relates the angular diameter distance to the last scattering surface with the 
angular scale of the first acoustic peak in the WMAP7 power spectrum, and is given by (for k = 0) (3lj 

1090 

R(n dni0 ,n b0 ,b) = [ dZ = 1.725 ±0.018, (12) 

J *H"dm0) "60j 0) 


where E(Cl dm o, ^60> — H(Q, dm o,Q b Q,b)/ Hq. It is important to point out that the measured value of R is model 
independent and, therefore, can be used as a test of cosmology. In general, the shift parameter alone is insufficient 
for this purpose. In ref. [32j it is shown how to derive the R parameter in a model independent way. In order to 
include the CMB shift parameter into the analysis, it is needed to integrate up to the matter-radiation decoupling 
(z ~ 1090), so that radiation is no longer negligible and it should be properly taken into account. For our model we 
used the minimization of the \ 2 statistic according to 

2 _ [1.725 - R(n dm0 ,n b0 ,b)} 2 

Xr - ooiP • (U) 

We used il b o = 0.04. In the Fig. 7. one can see the results for our model. It is clear that this test alone can not 
severely restrict the space of parameters. In particular, this test does not allows us to determine an upper limit for 
the parameter b for fixed O's. 

B. BAO 

The primordial baryon-photon acoustic oscillations leave a signature in the correlation function as observed by 
Eisenstein et al. [33| . This signature provides us with a standard ruler which can be used to constrain the following 
quantity 



A(Qdm0)^b0>b) — yQmJQ \_E(Z\, ridm0> ^60, &)] 
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where the observed value of A is A obs — 0.469 ± 0.017 and z\ = 0.35 is the redshift of the SDSS sample. The x 2 
statistics may be computed from 

a , q n [0-469 -A(n dm0 , 0.044, b)f 

XBAo( il dmO,lho,b) = 0017 1 ' ^ ' 

Here we will use the Q b o = 0.044 value to build the regions confidence. The results are presented in Fig. 8. We can 
see that BAO does not impose an upper limit for the parameter b, but it restricts the value of fidmo- 

C. H(z) 

The observational Hubble parameter depends on the redshift z in the form 

H ^ = -YTz% ^ 

Then, H(z) can be obtained as far as dz/dt is known. Jimenez et al. [34[ demonstrated the feasibility of the method 
while Simon et al. [35[ and Stern et al. [36[ obtained H(z) in the range of 0.09 < z < 1.8. Another way to get values 
of H(z) can be using the BAO peak position as a standard ruler in the radial direction. Using this prescription in 
the reference J33 the authors determined H(z) in z = 0.43 and z = 0.24. We employ the eleven data of [35ll36j and 
two data of [371 ] . The best fit values for the model parameters from observational Hubble data are determined by 
minimizing the quantity 
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Here we use the value flbo = 0.04. The free parameters are fidmO: Ho an d b. However, in our case, we have found that 
the one-dimensional probability distribution function in H is quite narrow and the values H are well approximated 
by the value Ho = 70.5 (minimizing the Xh considering the three free parameters). Therefore, one has only two free 
parameters, namely b and and the f2dm0 parameter. The results for this test are presented in Fig. 9. In general, this 
test provides the results which are similar to the ones of BAO. In particular, this is because the parameter b remains 
weakly constrained, but the restrictions are strong for the value of fWiO- 



D. Compilation Union2 



The supernovae la data give us the distance modulus \i to each supernova, that is given by 

r d L 



ji = m — M = 5 log 



25, (18) 



-Mpc\ 

where M is the absolute magnitude. The distance modulus can also be written as 

H = 5log 10 D L (z) +fio , (19) 

where Dl — Ha ^ L is the Hubble-free luminosity distance and is the zero point offset (which is an additional 
model-independent parameter) defined by 

Mo = 5 log 10 (^-) + 25 = 42.38 - 5 log 10 h . (20) 

In the present paper we used the Union2 dataset including 557 data of Amanullah et. al (38|, that includes the 
intermediate-z data observed during the first season of the Sloan Digital Sky Survey (SDSS)-II supernova survey [39| 
and the high z data from the Union compilation (4pj . 

In our case the XsNia is given by 

2 , s_ V- {iHheiPhZi) ~ Hobsizi)) 2 . . 

XSNIaKPi) — 2 ■ \ Zl > 

i—1 ® obs,i 

The x 2 function can be minimized with respect to the /j,q parameter, as it is independent of the data points and the 
dataset. The a b s are provided in j38| and include both the observational and intrinsic magnitude scatter. We assumed 
a Gaussian and uncorrelated probability distribution for the data. Expanding the equation above with respect to 
we obtain: 

X 2 (p l )sNia = A{ Pi ) - 2fx B( Pi ) + n 2 C( Pt ) , (22) 
which has a minimum for — B(pi)/C(pi), giving 

- XsNIa - MPi) - C{p~i) ' *~ 23 ^ 

where 

AOpA = Y^th-HoOs(Pi,»o = 0)) 2 (24) 

% 

B(pi) = , (25) 

i 

Cipi) = 4- (26) 
°i 

Let us note that the new Xsnio, i s independent on and can be minimized with respect to the parameters of our 
theoretical model. The results are shown in the Fig. 10. We can see that the parameters b and f^dmo are relatively 
constrained. In particular, the parameter b has a maximum likelihood value for b w 0.4 (see Table 1). 



8 



E. 2dfGRS 

We used the 2dFGRS data [4l[ to compare the power spectrum of our model and constrain the free parameters. 
To do this, we minimize % 2 defined as 



X2dFV l dmO, "60 > »J - 



07 



(27) 



Here P(fc) is given by equation ([SJ and the number of free parameters is two, since we fixed the value of fl^o = 0.04. 
In Fig. 3 we present the PDF for our model. The impact that has the model on linear power spectrum is dominated 
by the b parameter that is associated with the speed of dark matter particles. In the left figure we can see that if 
O(j m o is close to ACDM value, then, the b value is close to zero. But, if the value of b is, for example, 10 -4 , then the 
dark matter fraction decreases significantly. The corresponding data can be seen in Table 1. 




FIG. 3: The probability density function using 2dFGRS for Q m o (right) and 6(left). The confidence level with 1-a, 2-a and 3-cr. 
To build the PDF we marginalize the free parameters considering the intervals: fidmo e [0.05, 0.95] ACDM and b e [0.001, 0.4], 
for each case. As we can see the model is included, i.e, the first figure shows that for a fidmo ~ 0.25 corresponds to a b ~ 0. 



F. Combining the Datasets 

We considered that the observational data are independent, so we defined the Xtotai as 

Xtotal = XSNIa + xIaO + x\ + Xh + xldF ( 28 ) 
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TABLE I: In this table we compile the results of observational constraints using different dataset. The statistical error is 
68.3%(1ct) and the results in brackets correspond to an error of 2a. 



Data 


Xmin 


SldmO 


6 


2dFGRS 


21.495 


n nnfiq +(, ' 1604[+0 ' 34001 

U-UUU -0.0069[-0. 00691 


U ' UU -0.00l[-0.0011 


SNIa 


542.45 


n 2i + 002 


0A8t° o f o 


R+BAO+H(z) 


8.449 


0.28l°;gJ 


^ n +0.006 
~ u -0.000 


R+BAO+H(z)+SNIa(Base) 


547.383 


0.26t° :» 


U.UUUO±_ () 00064 


R+BAO+H(z)+2dfGRS 


24.24 


0.28l° :S 


^ n +0.001 
~ u -0.000 


R+BAO+H(z)+SNIa+2dfGRS 563.271 


o.26±g : gS 


~ u.uuui_ 0001 



The best fit values the model can be determined by minimizing the total \ 2 . For Gaussian distributed measurements, 
the x 2 function is directly related to the maximum likelihood estimator. The likelihood function is determined as 

L({l dm0 , b ) - L cxp ^-^LsL^j , (29) 

where Lq is a normalization constant. In order to constraint the parameters of our interest, we marginalize over the 
other parameters. In the Figs 4, 5 and 6. We displayed the results the \ 2 total. An important issue that needs to be 
considered is the difference between the PDFs of the background tests (BAO+R+H+SNIa) and linear perturbations. 
In the case of fidmo parameter, we see that the low probability in the range of fidmo < 0.1 eliminates high probability 
of this sector in case of perturbations (Fig. 3). Therefore, the result are close to ACDM model. Effects of structure 
formation may be critical to study the characteristics of models that include warm dark matter, where the speed of 
the particles is no longer negligible as in the case of CDM. In the table 1 we compile the results of observational 
constraints at la and for the power spectrum, we give also the estimations at 2a. 
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FIG. 4: Confidence regions at 1-cr, 2-a and 3-<r levels from inner to outer respectively on the (fidmo , b) plane for our relativistic 
model in the flat case. The effect of type la supernovae is not apparent when we consider only tests of background (H+R+BAO), 
but it is important when combined with power spectrum data, for example, the constraints on fidmo parameter are more narrow 
and the upper limit for b parameter decreases significantly, an order of magnitude, compare the last two figures. 



RESULTS AND DISCUSSION 



The Reduced Relativistic Gas model (RRG) enables one to determine the equation of state which interpolates 
between the radiation-dominated and the matter-dominated eras. It can represent a warm dark matter, characterized 
by the parameter b, where v ~ be. The main difference from WDM and CDM come from the behaviour at small 
scales, even during the relativistic phase. WDM is composed from lighter particles with respect to the CDM particles, 
and density fluctuations are suppressed at small scales, presenting the effect of a free stream. In our case, such small 
scales effect are not relevant, since we are interested in scales corresponding to the linear regime, and transition is 
smooth, not taking into account the details of the decoupling of the warm dark matter particles. In such a case, the 
main effect of our model is due to the non-zero, even though small, value of the velocity of the particles. For our 
purposes, we can ignore these small scales process and take into account only the informations encoded in the transfer 
function for our relevant scales. For very large scales, the WDM and CDM transfer functions are essentially the same, 
and small differences appear only at intermediate scales, as we could expect. This is in agreement with our results 
which are essentially the same for the WDM and CDM transfer function. Hence, the dark component would have a 
small (but non zero) velocity today. 

The model has been tested using four background observational tests: Supernova type la (Union2 sample), H(z), 
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FIG. 5: The one-dimensional probability distribution function (PDF) for the OdmO parameter. The parameter b was marginal- 
ized in the range [0, 1]. 

CMB (R factor) and BAO. Moreover, a detailed study of structure formation at linear level has been performed using 
the 2dFGRS data for matter power spectrum. The different tests have been crossed in order to obtain a more clear 
prevision for the free parameters, which are essentially the velocity parameter for the dark matter particles b and the 
dark matter ratio to the critical density fidmo- All the analysis has been performed using the flat universe prior. 

In general, the background tests predict a f^dmo very similar to the value obtained in the ACDM model, except 
for the CMB R factor which predicts a universe essentially dominated by dark matter. For the velocity parameter b, 
the background tests predict a maximum probability for a finite value at b -C 1, except for BAO, where a large value 
for b is predicted. But the most striking feature of the observational tests concerns the formation of structure, for 
which the maximum probability for Sldmo occurs at a zero value. This seems to be a consequence of the restriction 
of the analysis to a linear level, since a certain amount of dark matter is necessary in order to have the formation 
of structure process. In any case, the results concern a probability distribution: a small amount of dark matter 
is certainly admitted, but much less than that predicted by the ACDM model. For the parameter b, the PDF is 
concentrated around zero, but there is another region of much smaller, albeit non-zero, probabilities from b ~ 0.04 
on. We remember that the estimations of warm dark matter in dwarfs spiral galaxies predict b ~ 10 -3 [42j]. In this 
context, our model seems to be a possible alternative for understanding the formation of structures. Additionally as 
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we had seen in the RRG the dispersion velocity is a free parameter predictable directly and not derived. 

The crossing of the background tests with the ones related to cosmic perturbations leads to scenarios similar to 
the ACDM model but with a quite large dispersion in the estimations for b. Here we must distinguish some special 
features. When the SNIa test is not considered, the range of allowed value for b becomes larger. The addition of the 
structure formation to the background tests does not change considerably the scenario, in spite of the fact that test 
alone would lead to a very small value of fidmo- This occurs because for small Sldmo the PDF is essentially zero for the 
background tests: the crossing of perturbative and background tests renders this possibility very unlikely. However, 
considering an error of 2a (see table 1) our value of QdmO includes the value of the ACDM model Sldmo = 0.229 ±0.015 

Em 

One can resume the main results of this work in the following way. The background tests result, without the SN 
data, is centered in the usual values characterizing the ACDM model, that means b w and fidmo ~ 0.26, but with a 
non-negligible dispersion. The inclusion of the SN data sharpens the confidence region, leading to stronger evidences 
in favor of the ACDM model. In both cases, the inclusion of the LSS test enlarges the dispersion mainly in the 
parameter b direction. For the one-dimensional estimations, those related to the dark matter density changes very 
little with, while for the b parameter, the inclusion of SN in the background tests compress the PDF distribution near 



13 



b = 0, and the inclusion of the LSS tests accentuates this effect. This seems to be in contradiction with the previous 
two-dimensional description, but it must be seen as an effect of the marginalization on the density parameter. 

In the literature the mechanisms most used for to study warm dark matter has been based in a distribution 
function of the form f x (v) = eP / a T y+1 43] where T 7 is the photon temperature, v — p/(p 2 + m 2 ) 1 / 2 , and three 
parameters a, /?, and m x . The power spectrum today and the relation between pressure and energy density is given 
with a function of this distribution. In this context the transfer function can be approximated by the fitting function 
T(k) = [1 + (afc) 2 - 24 ] -4 - 46 [29( where a is a function of free parameters of WDM. We have used this function with 
m x = IkeV and m x = 0.25keV and observe that the obtained results are very close to those obtained with the BBKS 
function, therefore we our results arc not affected significantly by the scale of mass. 

It is important to remember that the standard model, ACDM, appears to be in conflict with observations on 
subgalactic scales. There are two major conflicts between ACDM models and observations in the Local Universe. 
First, the inner mass density profiles of simulated dark matter halos are more cuspy than inferred from the rotation 
curves of dwarfs and low surface brightness galaxies. Currently, this conflict seems well established [HI]. Therefore, the 
ACDM models in its current form cannot provide the unknown mass that surrounds galaxies. The second problem 
is related with the N-body simulations of ACDM models that predict large numbers of low mass halos greatly in 
excess when compared the observed number of satellite galaxies in the Local group. These and other issues suggest 
that it is necessary to investigate different theoretical possibilities. In this paper we have shown that, in part of 
linear perturbations alone, our RRG model for a WDM provides results which are different from the standard ACDM 
model. It is possible, in principle, that the parameter estimations from other observational data can be modified by 
the quantum corrections in the form discussed in (24| and pBj , such that eventually we may converge to an alternative 
cosmic concordance model. The present investigation can be seen as a first step towards this program. 

It is clear that to better quantify these qualitative expectations, one needs to take into account the nonlinear effects, 
including the dynamics of the collapse and the halo concentrations, e.g., by using algorithm suggested in Ref. [45j . 
etc. Another important thing to do is to consider the other sets of observational tests, like the integrated Sachs- Wolfe 
effect, in the first order in quantum corrections and see whether they are compatible with the small flodm and small 
WDM model. 

Acknowledgment s . 

The work of J.F. has been supported in part by CNPq and by FAPES. The work of A.M.V.T. and I.Sh. has 
been supported in part by CNPq and FAPEMIG, I.Sh. also was supported by ICTP Senior Associated Membership 
Program. 



14 



Appendix 





FIG. 8: The PDFs using BAO. The confidence level with 1-a, 2-cr and 3-a(left) 




FIG. 10: The PDFs using SNIa Union 2. The confidence level with 1-cr, 2-cr and 3-a(left) 
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